Magnetic quantum oscillations in YBa 2 Cu 3 06.6i an d YBa 2 Cu 3 06.69 m fields of up to 
85 T; patching the hole in the roof of the superconducting dome 
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We measure magnetic quantum oscillations in the underdoped cuprates Yl^CusOe+a; with x — 
0.61, 0.69, using fields of up to 85 T. The quantum-oscillation frequencies and effective masses 
obtained suggest that the Fermi energy in the cuprates has a maximum at p ~ 0.11 — 0.12. On 
either side, the effective mass may diverge, possibly due to phase transitions associated with the 
T — limit of the metal-insulator crossover (low-p side), and the postulated topological transition 
from small to large Fermi surface close to optimal doping (high p side). 



PACS numbers: 71.18.+y,74.72.Gh,71.45.-d 

One of the significant landmarks in the study of the 
"High-T c " cuprates is the observation of Shubnikov-de 
Haas and de Haas-van Alphen oscillations [TJ HJ [3J 2J 
03 [5J [7] in high magnetic fields. Such magnetic quan- 
tum oscillations (MQOs) are the signature of a Fermi 
surface (FS), and their temperature (T) and field (B) 
dependence suggest a relatively conventional Fermi liq- 
uid [TJ [31 IH |3J 13 [5] , rendering some theories of the 
cuprate normal state untenable [SJ. Though there are 
attempts to explain the MQOs using more exotic mod- 
els pUl [TT| [12] , these seem unable to describe aspects of 
the data (e.g. multiple MQO frequencies, MQOs periodic 
in 1/B, realistic effective masses). 

However, published MQOs cover only a restricted re- 
gion of hole doping p. In particular, data on under- 
doped YBa 2 Cu306+£c correspond to 0.49 < x < 0.54 
(0.0925 < p < 0.10) [TJ EJ. As this is also the x 
range blighted by the ortho-I/ortho-II structural insta- 
bility [T5], it is natural ask whether the observed FSs are 
a consequence of, or related to, this phase separation. 
Moreover, the only higher-p data for the underdoped 
side of the superconducting dome are for YBa2Cu40g 
(p w 0.125 — 0.14 [3j[4]). These may be untypical because 
of the different crystal structure. Here we therefore re- 
port MQOs in the underdoped cuprates YBa2Cu306.ei 
(p w 0.11) and YBa 2 Cu 3 6 .69 (p » 0.125). We find 
that both exhibit a dominant MQO frequency F » 
550 — 570 T, similar to the a frequency observed in the 
0.49 < x < 0.54 samples [TJ [5J [5J [7J . On close examina- 
tion, the p ss 0.11 sample exhibits additional MQO fre- 
quencies, some attributable to warping of the FS due to 
a finite interlayer transfer integral. Effective masses m* 
found for both compositions are less than 2.0m c , lighter 
than their equivalent in YBa2Cu40g (SJH]. 

Single crystals of YBa2Cu30g.6i and YBa2Cu30g.69 
are grown and oxygenated as described before [TJj . Sam- 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




as 

10 20 30 40 50 60 70 80 

Magnetic Field, B, (T) 



FIG. 1: (color online.) PDC frequency / versus field B 
for single crystals of YBa2Cu3 06.69 (upper blue trace) and 
YBa 2 Cu 3 6 .6i (lower red trace); T = 1.5 K for both. The 
drop in / corresponds to the irreversibility field [U [5] . MQOs 
are clearly visible at high fields for YBa2Cu306.6i- The inset 
shows the B versus time t profiles for the 60 T Long-pulse 
and 85 T Multishot magnets. 

pies are polished to sizes 0.3 x 0.3 x 1.5 mm 3 , with the 
long axis parallel to c. Compositions arc inferred by mea- 
suring T c in a SQUID magnetometer, and using the p and 
x versus T c relationships given in Ref. [13]. The MQO 
experiments employ the same system as in Ref. [J ; a coil 
of 8-15 turns of 44 or 50-gauge Cu wire is wound around 
the sample, the planes of the turns roughly perpendicu- 
lar to c. The coil is part of a tank circuit driven by ei- 
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ther a tunnel-diode oscillator (TDO) 15J or a proximity- 
detector circuit (PDC) [16 ; shifts in resonant frequency 
/ are caused by changes in the skin-depth (normal state) 
or penetration depth (superconducting state) [15] ■ No 
significant differences are noted between PDC and TDO 
data. A heterodyne system measures /; the oscillator 
output is mixed down using two mixer/filter stages to 
about 1 MHz and the resulting signal digitized directly at 
10 7 samples/s using a National Instruments PXIe-1062Q 
digitizer. Fields are provided by the 85 T Multi-shot 
(MSM) and 60 T Long-pulse magnets at NHMFL Los 
Alamos |H H] and a 65 T short-pulse magnet at Oxford. 
The purpose of the range of dB/dt (~ 100 - 15000 Ts" 1 ) 
is to characterize and eliminate the effects of sample 
heating due to induced currents and dissipative vortex 
motion [4 . The field is measured using a pick-up coil 
calibrated by the belly MQOs of Cu [T7]. Four crys- 
tals of YBa2Cu306.6i and two crystals of YBa2Cu306.69 
are studied; results are consistent between crystals of the 
same x and between different magnets. 

Fig. [I] shows data for YBa 2 Cu306.6i and 
YBa 2 Cu 3 6 . 6 9 measured in the 85 T MSM at T = 1.5 K; 
samples are heat-sunk to a sapphire chip and immersed 
in 4 Hc liquid. The frequencies are obtained by Fourier- 
transforming the signal using a moving time-window 
20 /is long, and then adding the offset removed by the 
mixers. The prominent drop in / around 25 T [x = 0.61) 
or 35 T (x = 0.69) is attributed to the irreversibility 
field |H |5]. Above this, features are discerned in the 
data, corresponding to Shubnikov-de Haas oscillations 
in the conductivity [U [5] . Owing to the proportionality 
between change in conductivity and shift in / |15j . the 
conductivity MQOs give oscillations in /. 

We first turn to YBa2Cu 3 06.6i for which MQOs are 
visible in the raw data (Fig. [I]); below, we see that MQOs 
are less prominent in YBa2Cu30g.69 due to a shorter ap- 
parent scattering time, r. To make the MQOs more visi- 
ble, the slowly- varying background due to the semiclassi- 
cal magnetoresistance is removed by subtracting a third- 
order polynomial in B. Figure [2]ja) shows some resulting 
MQOs for YBa2Cu306.6i! here, random noise from the 
power supply of the 85 T MSM [5] is mitigated by av- 
eraging three upsweeps and three downsweeps and then 
smoothing using a Savitsky-Golay routine. The result- 
ing data exhibit MQOs above about 40 T. On Fourier- 
transformation, the dominant peak is at F = 570 T 
(Fig. (2jb)), similar to the so-called a MQO frequency 
(500 - 550 T) in 0.49 < x < 0.54 samples [I] IB] and 
the dominant frequency in YBa2Cu40s OHj. 

However, in the case of a single extremal FS cross- 
section, one would expect the MQO amplitude to grow 
uniformly with increasing field [19]. The MQOs in 
Fig. [2ja) do not do this; they are modulated by what 
appears to be a beat frequency, a phenomenon noted 
other cuprates [5] [3] [TB] . This is also seen in the Fourier 
transform (Fig. [2^b)), where the peak at 570 T is obvi- 
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FIG. 2: (color online) (a) PDC resonant frequency / for a 
YBa2Cu3Cv6i crystal after background subtraction to leave 
the oscillatory component A/; the trace (thick black curve) is 
an average of three magnet sweeps (T = 1.5 K). The thinner 
red line is a fit of Eq.[l]for MQO frequencies 589 T and 479 T 
with r ss 0.07 ps. (b) Fourier transform of the data in (a) 
(black thick curve); the large peak is centred on 570 T. The 
thinner red curve is the sum of two gaussians (fine green lines) 
at 466 and 593 T. (c) Plot of MQO amplitude A divided by 
T versus T; diamonds are from the upsweep of the field and 
dots from the downsweep. The curve is a fit of Eq. [2] giving 
m* = 1.6 ± 0.1m c . (d) Residual (i.e. (data)-(fit)) from (a) 
versus field (black thicker line) . The thinner red curve is a fit 
of Eq. [l]for a single MQO with F = 270 ± 20 T. 

ously asymmetric and may be fitted by two overlapping 
Gaussians centred on 466 ± 10 and 593 ± 5 T. The pres- 
ence of two relatively closely-spaced Shubnikov-de Haas 
frequencies with similar amplitudes is suggestive of the 
beats caused by "neck and belly" oscillations of a quasi- 
two-dimensional FS that is warped due to a finite inter- 
layer transfer integral t£ 0HH]- To model this, we sum 
two components of the Lifshitz-Kosevich formula [HI [T3] 
with MQO frequencies F\ and F 2 , amplitudes a\ and a 2 
and phases <pi and <p2 '■ 
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Here m* is the effective mass, and r _1 is an effective 
scattering rate; we assume that m* and t are the same for 
the neck and belly oscillations. Independently, m* may 
be constrained by the way in which the amplitude of an 
individual MQO, or the Fourier amplitude of a transform 
over a restricted field range varies with T; 
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where A is the amplitude. All of the fits (e.g. Fig.[2](c)) of 
individual MQOs or Fourier amplitudes yielded m* val- 
ues in the range 1.5 — 1.7m e , irrespective of sample, field 
range or sweep rate, leading us to m* — 1.6±0.1m G . Hav- 
ing constrained m* , a fit of Eq. [l]to the data (Fig. [2j») 
yields MQO frequencies 589 ± 5 T and 479 ± 5 T, and 
r « 0.07 ps. These values are close to those obtained in 
the two-gaussian fit of the transform in Fig. [2^b) . 

Beside the peak at 570 T in the transform (Fig. [2](b)), 
there is a feature at F ss 250 T. This appears to corre- 
spond to an actual MQO series, as is seen by subtracting 
the fitted Eq. [I] in Fig. [2^ a) from the data. The resid- 
ual (Fig. |2jd)) is oscillatory, with a direct fit yielding 
F = 270 ± 20 T, close to the value suggested by the peak 
in the Fourier transform [5D]. Unfortunately, the MQOs 
are too poorly defined to permit estimates of m* or r. 

To summarize for YBa2Cu306.6ii our data suggest 
three FS cross-sections, with MQO frequencies 270 T, 
479 T and 589 T; other peaks in the transform at higher 
frequencies (Fig. |2^b)) are attributable to harmonics of 
these [H]. The 479 T and 589 T MQOs are likely 
the neck and belly oscillations of a warped quasi-two- 
dimensional FS, with m* = 1.6 ± 0.1m c ; this is probably 
the equivalent of the dominant a frequency in other un- 
derdoped cuprates [U |2j [5j [5J . The frequency difference, 
AF w 110 T, between neck and belly oscillations sug- 
gests [IS] an average interlayer transfer integral t^r = 
hAF/(4m*) = 2.0 ± 0.1 meV for YBa 2 Cu 3 6 . 6 i, higher 
than the values 1.4—1.7 [23] obtained for YBa2Cu306+ K 
(x = 0.51, 0.54) [2j [18]. This increase in with p is not 
unexpected; the lattice parameter c declines with p [13] . 

Figure [3^a) shows an example of the MQOs observed 
in YBa2Cu 3 06.69- In contrast to YBa 2 Cu306.6i, where 
MQOs appear around 40 T (Fig. (2j»), the MQOs 
here are not distinguishable from the background un- 
til about 60 T (Fig. ^a)). The non-sinusoidal appear- 
ance of the MQOs again suggests the presence of more 
than one frequency, but sadly, the limited field window 
over which MQOs are seen both precludes a "neck and 
belly" analysis (Eq. [TJ Fig. [2^a)) and limits the resolu- 
tion of a Fourier transform. Instead, we plot MQO in- 
dex versus 1/2? in Fig. [3|b) to find a mean frequency 
of 550 ± 20 T [24 . Fitting the MQO amplitudes versus 
T for YBa2Cu306.69 yields m* — 1.8 ± 0.3to c , similar 
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FIG. 3: (color online) (a) PDC resonant frequency for a 
YBa2Cu30e.69 crystal after background subtraction to leave 
the oscillatory component A/; the trace is a smoothed average 
of three magnet sweeps (T = 1.5 K). (b) Oscillation index 
versus reciprocal magnetic field for the MQOs in (a) (points) ; 
dips in A/ are indexed by integers and peaks by half integers. 
The straight line is a fit with a gradient of 550 ± 20 T. 
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FIG. 4: (color online) (a) Summary of MQO frequencies 
versus p for underdoped cuprates: for YB^Cu^Oq+x, < are 
from Ref. |5j, from Ref. [2], and o from this work; T c and 
x values are converted to p using Ref. [13]. YBa2Cu408 data 
from Refs. [3 ( 4 are squares; the horizontal bar is the spread 
in p values given for YBa2Cu40s [3] [4]. Solid symbols (e.g. 
• ) show the dominant (F a ) frequency obtained from Fourier 
analysis; open symbols are from more detailed analysis (e.g. 
Figs. [5Fa,d) or Refs. [2] [18]). (b) Effective mass of the domi- 
nant (a) MQO frequency F a versus p\ symbols are the same 
as in (a) except are from Ref. [?]• (c) Fermi energy from 
F a and m*; for ps where several values are given, we take the 
average. Points are data and the curve is a parabolic fit. 

to the 1.6 ± 0.1m o for the analogous MQO frequency in 
the YBa 2 Cu30 6 . 6 i (Fig. [2^c)). A Dingle analysis (i.e. 
a plot of log e (AB5 sinh(14.69m*T/23) versus 1/2? [19], 
where A is the oscillation amplitude) yields r 0.04 ps, 
~ 2 times smaller than that for YBa2Cu306.6i- This ac- 
counts for the higher fields required to observe MQOs 
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in YBa2Cu3C>6.69- Ref. ES] attributes the dominant 
Landau-level broadening to quasistatic spin disorder also 
observed in neutron experiments and parameterized by 
a correlation length £ [57] HE]- It is notable that £ 
decreases with increasing p H7J [55], and this may 
account for the shorter r of the x — 0.69 samples. 

Figure [4] compares the data obtained here with simi- 
lar results from other underdoped cuprates, all of which 
have a dominant MQO frequency F a ps 500 — 660 T. 
Figure |4ja) shows both F a and other MQO frequen- 
cies < 1000 T that have been resolved (this paper, 
Refs. [2J HSJ). If we attribute F = 540 and 450 T for 
p ps 0.10 [2] and F = 590 and 480 T for p = 0.11 (this 
work) to the belly and neck oscillations of the a Fermi 
pocket, then there seems to be a trend, continued by 
YBa2Cu308, for the a pocket to grow with rising p. It 
also seems that samples from the orthoI-II region are un- 
exceptional, continuing the trend seen in this work to 
lower p. The weaker MQOs with F = 630 T (p = 0.975, 
0.10) [2 F — 270 T (p = 0.11) are qualitatively similar 
to extra pockets predicted by FS reconstruction due to 
various types of symmetry breaking; e.g. an incommen- 
surate spin-density wave [29] produces a plethora of FS 
sheets, both smaller and larger than the a pocket, whilst 
a pocket with F w 250 T is an explicit prediction of in- 
commensurate d-density-wave order [30, 31]. Meanwhile, 
the a effective masses show a "bowl" -shaped dependence 
on p, with a minimum at p ps 0.11. 

To vizualize effect that these changes have on the car- 
rier system, Fig. Etc) plots the effective Fermi energy 
Ep for the a pocket, Ep = hF a /m*, using data from 
Figs. EJa,b). It seems that the Fermi energy reaches 
a maximum at p sa 0.115, but decreases either side of 
this, suggesting that m* may diverge at p p=> 0.087 and 
p rj 0.14, the latter p being poorly constrained by the ex- 
isting data [33J . The lower p value suggests the point at 
which Metal-Insulator transition tends to T = [8l 132] . 
The upper may signal the topological transition from 
small to large FS thought to occur close to optimal dop- 
ing [U [BJ, though experimental confirmation of an un- 
reconstructed FS in overdoped YBa2Cu306+a: is as yet 
lacking. By analogy with heavy-fermion superconduc- 
tors 134] . both of the m* divergences may represent 
quantum critical phase transitions. 

In summary, we report MQO frequencies and effective 
masses for the underdoped cuprates YBa2Cu30g+a; with 
x = 0.61, 0.69, filling in a considerable gap in the FS 
versus hole doping diagram. In conjunction with other 
data, our results suggest that the Fermi energy may reach 
a maximum around p sa 0.11 — 0.12, and collapse on ei- 
ther side due to divergence of the effective mass. The 
divergences are perhaps associated with quantum-critical 
phase transitions associated with the T = limit of the 
metal-insulator transition (low-p side), and the topolog- 
ical transition from small to large FS close to optimal 
doping (high p side) . 
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